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ABSTRACT 

We report new deep ALMA observations aimed at investigating the [CII] 158yum line and con¬ 
tinuum emission in three spectroscopically confirmed Lyman Break Galaxies at 6.8<z<7.1, 
i.e. well within the re-ionization epoch. With Star Formation Rates of SFR ~ 5 - 15 M© yr“^ 
these systems are much more representative of the high-z galaxy population than other sys¬ 
tems targeted in the past by millimeter observations. For the galaxy with the deepest obser¬ 
vation we detect [CII] emission at redshift z=7.107, fully consistent with the Lya redshift, 
but spatially offset by 0.7" (4 kpc) from the optical emission. At the location of the optical 
emission, tracing both the Lya line and the far-UV continuum, no [CII] emission is detected 
in any of the three galaxies, with 3cr upper limits significantly lower than the [CII] emission 
observed in lower reshift galaxies. These results suggest that molecular clouds in the central 
parts of primordial galaxies are rapidly disrupted by stellar feedback. As a result, [CII] emis¬ 
sion mostly arises from more external accreting/satellite clumps of neutral gas. These findings 
are in agreement with recent models of galaxy formation. Thermal far-infrared continuum is 
not detected in any of the three galaxies. However, the upper limits on the infrared-to-UV 
emission ratio do not exceed those derived in metal- and dust-poor galaxies. 
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1 INTRODUCTION 


Millimeter and submillimeter observations of distant galaxies are 
a powerful tool to trace the evolution of the gas and dust content 


(e.g. Tacconi et al.|2013| 

Santini et al. 

2014 

Bothwell et al.|2013[ 

IGenzel et al.|201ol|Dadd 

li et al.|2010 

|Cari 

& Walter|2013J, to 


investigate their dynamics ( e.g.|Wang et al.|2013[|De Breuck et af] 


20T4l|Carniani et al.|2013||Carilli et al.|2013||Willott et al.|2013 

2015a|l to trace obscured star formation (e.g.|Swinbank et al.|2014 

Gallerani et al.||2012| 

INegrello et al.||2010| |Hezaveh et al.||2013 

Schaerer et al.|l2015a| 

►, and even to trace AGNs and their effect 

onto their host galaxy 
IGallerani et al.|20141). 

(|Maiolino et al.||2012| |Cicone et al.||2015| 


Molecular (mostly CO) transitions are used to trace the con¬ 
tent of molecular gas, while continuum emission is used to trace 
warm dust heated by ongoing star formation (and possibly by an 


AGN). At z>4 mm/submm observations are even more effective 
since the luminous [CII]158/im far-infrared fine structure line is 
redshifted into the primary bands of atmospheric transmission, at 
frequencies v < 350GHz (|Maiolino et al.|2005||2009||Walter et al.| 
|205^|De Breuck et al.|2014 . The [CII]158/<m line is one of the 

main coolants of the ISM, both in local and high-z galaxies, gen¬ 
erally excited by heating of UV photons in star forming regions, 
primarily associated with their Photo-Dissociation-Regions (|Mad- 


Mil 


den et al.||1997[ Kaufman et al.||1999[ |Gracia-Carpio et al.||201 


Pineda et al.|2014| l, although contribution from diffuse atomic and 
partially ionized gas, as well as from shocked gas, has been in¬ 
ferred in some galaxies ( Pineda et al.|2014[ Appleton et al. Pol'S! 


[Velusamy & Langer|2014[ Cormier et al.|2012|l. 


Until recently, at z>4 [CII], CO and continuum emission 
had been detected only in extreme galaxies such as quasar hosts 




































2 Maiolino et al. 


Table 1. Sample and parameters of resulting from the observations. 


Name 

RA(J2000) 

[deg] 

DEC(J200) 

[deg] 

^Lyor 

SFR" 

[Moyr-i] 

yobs([CII])'’ 

[GHz] 

beam^ 

[min"xmaj" (PA°)] 

U(int.)^ 

[hours] 

antennae 

^ cont 

I/'Jy] 

^line 

[mJy] 

L([CII])« 

[10^ Lo] 

BDF-3299 

337.0511 

-35.1665 

7.109'" 

5.7 

234.374 

0.55 X 0.76 (85) 

5.1 

25-36 

7.8 

0.062 

<2.0 

BDF-512 

336.9444 

-35.1188 

7.008'" 

6.0 

237.330 

0.72 X 0.50 (68) 

1.4 

29 

17.4 

0.171 

<6.0 

SDF-46975 

200.9292 

-^27.3414 

6.844"' 

15.4 

242.290 

0.96 X 0.82 (79) 

2.0 

40 

19.2 

0.173 

<5.7 


Notesi'^This is the SFR inferred from the rest-frame UV continuum and is adapted to the new calibration given in |Kennicutt & Evans] j20 12) , who adopt a 
Kroupa IMF. ^Expected redshifted frequency of [CII] (LSRK) from the Lya redshift. ‘^This column gives the synthesized beam (median value across SWl) 
minor and major axes (in arcsec) and position angle (in degrees, anticlockwise relative tot the North). ^On-source integration time (in hours). ^Rms on the 
continuum. ^Rms measured in the spectrum extracted from the central beam, avoiding regions of higher noise associated with atmospheric absorption (except 
for SDF-46975, see text), in spectral channels of 100 km /s. ■^The upper limits o n the [ CII] luminosities a re at 3cr, and are calculated on a width of 100 km/s. 

^From f^nzella et aL]j201l) . ^From |Ono et al.| ( |201^ . 


or SMGs, characterized by star formation rates of the order of 
1000 Mq yr"\ certainly not representative of the bulk of the galaxy 
population. In the last few years a few galaxies at z~4-6 with more 
modest star formation rates (20 - 300 M© yr"^) have been detected 
in [CII] or in continuum ( Carilli et al.|2013[[W agg et al. 2012[|Car-| 


|niani et aT||2013[ | Williams et al.||2014[ [Riechers et ai.||2014| ). It 


should be noted that even these SFR are not yet really representa¬ 
tive of the bulk of galaxyies in the early Universe, which have SFR 


typically lower than 10 M© yr ^ ( 

Salvaterra et al.[2011 

Finkel- 

[stein et al.[2012|[Dayal et al.|2014| 

Robertson et al.|2015[). 

At z~7. 


well within the re-ionization epoch, only one quasar host galaxy 
has been detected in [CII] and in continuum (at z=7.08, |Veneman^ 
|et aL||2012| ). “Normal” galaxies at z~7, with star formation rates 
more representative of the galaxy population at these early epochs, 
have not yet been detected either in [CII] or in continuum ( |Ota et al.| 
|2014[[Schaerer et al.|2015b||Gonzalez-L6pez et aL|2014| ). We also 

note that z~7.7 is the maximum redshift at which galaxies are spec¬ 
troscopically confirmed ( jOesch et al.|2015j ) and we recall that, even 
with the broad ALMA band, a spectroscopic redshift is needed to 
observe far-IR transitions redshifted into the millimeter band. 

We have undertaken an ALMA programme targeting a sample 
of three spectroscopically confirmed galaxies at 6.8 < z < 7.1 
with the aim of detecting or constraining both the [CII] line and 
the continuum thermal emission. 

For consistency with other papers on similar sources 
we assume the following cosmological parameters: Hq = 
70 km s"^ Mpc"^ Qa = 0.7 and = 0.3. 


2 SAMPLE SELECTION, OBSERVATIONS AND DATA 
ANALYSIS 

The list of galaxies observed is given in Tab. 1. These targets 
were initially selected as Lyman Break Galaxies through their z-Y 
dropout ( [Castellano et al.|2010[ [Ono et al.|2012[ l and all of them 
have solid spectroscopic confirmations, through the detection of 
Lyof with clearly asymmetric profile due to absorption of the blue 
side by the intervening IGM ( [Vanzella et al.[201 l1[Ono et al.[2012[ |. 

The redshift given in Tab. 1 refers to the redshift inferred from 
the peak of the Lycr emission. This is probably significantly in ex¬ 
cess (even by a few 100 km/s) with respect to the real systemic 
redshift of the sources, because of IGM absorption of the blue side 
of Lyof, and possibly also because of the effect of outflows (e.g. [Pet-[ 
[tini et al.|2002[ |. The quoted star formation rates are based on their 
UV emission, assuming no dust extinction (which however is likely 
to be very low given the very steep UV continuum, [Vanzella et al.[ 
[2011[ ). We have also corrected them, relative to the original papers. 


by adopting the new relation provided by [Kennicutt & Evan^ ( [2012[ ) 
(who use a Kroupa IMF). 

The ALMA observations were obtained between October 

2013 and April 2014, with a number of antennae ranging from 26 
to 40, depending on the specific observation (by also taking into ac¬ 
count the fact that a few noisy antennae were removed -flagged as 
“bad”- from further reduction). An early observation of BDF3299 
in May 2013 was not considered, because it was much noisier than 
the other data. The antennae were distributed in a semi-compact 
configuration, with most of them within about 200m from the array 
center, but with one or two antennae located as far as 500-1000m. 
The Precipitable Water Vapor during the observations ranged be¬ 
tween 0.9 mm and 1.8mm, with an average of about 1.2mm. 

Observations were performed in Frequency Division Mode 
(FDM). Out of the four spectral windows, SPWl was always cen¬ 
tered on the expected frequency of the [CII] line, in the Upper 
Side Band. This spectral band was set to a spectral resolution of 
2.5 km/s. SPWO was located on the continuum next to SPWl (on 
the higher frequency side), while SPW2 and SPW3 were located in 
the Lower Side Band to sample the continuum. The spectral reso¬ 
lution of the “continuum” spectral bands was different depending 
on the specific observation. 

The phases were centered at the nominal optical position of 
the source (but see discussion below about astrometric offset be¬ 
tween optical and millimeter data for BDF3299). 

The band pass calibrators for BDF3299 and BDF512 were 
J2056-4714 and J2258-2758, while for SDF46975 the bandpass 
calibrator was J1331-^3030. The flux calibrator for BDF3299 and 
BDF521 was Neptune, while the flux calibrator for SDF46975 were 
Ceres and Pallas. The phase calibrators for BDF3299 and BDF521 
were J2223-3137 and J2247-3657, while the phase calibrator for 
SDF46975 was J1058-h0133. 

We note that the observations of BDF3299 were obtained in 
two different runs, one in October/November 2013 and one in April 
2014. For the April 2014 observations, we moved the central fre¬ 
quency of SPWl by about 100 km/s with respect to the observation 
in October 2013 and also the phase center was offset by about 0.6". 
Therefore, before combining the data, the data obtained in April 

2014 were phase-shifted to match the phase center of the observa¬ 
tions in October 2013. 

The ALMA data were reduced and analyzed by using CASA 
version 4.2.1. Unfortunately this CASA version does not assign 
physically-based weights to data taken at different epochs, with dif¬ 
ferent integrations times, different atmospheric opacities and with 
different spectral resolutions. As a consequence, before combin¬ 
ing the visibilities of two observations at two different epochs, we 
manually re-scaled their relative weights so that they matched the 
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Figure 1. Top: spectra of the three sources in our sample, centered at the location of the optical counterpart, extracted from the central beam. The velocity 
reference is set to the redshift defined by the Lya peak. Bottom: noise in the same velocity range measured across each plane associated with each spectral bin. 
The regions where the noise is higher are associated with spectral regions of telluric absorption. The spectral binning of these spectra is 40 km/s. The vertical 
red dashed line in the leftmost panel indicates the velocity of the [CII] line detected next to BDF3299 in clump “A”, illustrating that the line is located in a 
region of low noise. 


noise variance measured in the two cubes. For the same reason, the 
same process had to be applied when combining data from differ¬ 
ent spectral windows (e.g. to generate the continuum maps) taken 
with different spectral resolution, i.e. the weights of the different 
spectral windows (when take with different resolutions) had to be 
adjusted manually. 

The cubes were created with natural weighting and by shift¬ 
ing the channels to a common velocity reference (i.e. by setting 
mode=’velocity’ in the CASA clean command). 

The resulting cubes images typically have a beam size of about 
0.7", yet the detailed beam shape and orientation for each obser¬ 
vation is given in Tab[^ No cleaning was performed in any cube, 
since there are no sources bright enough in the field of view. The 
only exception is the continuum image of BDF3299, in which a 
400yuJy serendipitous source is detected and which has enabled us 
to perform the cleaning (with 500 iterations) of the continuum map. 

The primary beam Full Width Half Maximum at this fre¬ 
quency is 22". 

The continuum maps were obtained by collapsing the four 
spectral windows (by taking care of the proper weighting be¬ 
tween spectral windows with different channel widths, as discussed 
above). 

The final sensitivity reached in each set of data is given in 
Tab.l, both for the continuum and for spectral bins (as detailed in 
the following). In the deepest cube (the one on BDF3299) we reach 
a continuum sensitivity of 7.8yuJy by using the equivalent full spec¬ 
tral band of 7.8 GHz. The sensitivity to line detection depends on 
the specific frequency, since each band is affected by atmospheric 
absorption in some spectral regions. The sensitivity also depends 
on the spatial location, both because of the shape of the primary 
beam, which decreases the sensitivity at large radii from the phase 
center, and because continuum serendipitous sources may gener¬ 
ate sidelobes in their vicinity. In Tab. 1 we give the rms in spectral 
bins of 100 km/s obtained by extracting the spectrum on the central 
beam (phase center) and taking the rms spectrally, in the spectral 
region close to the expected frequency of the [CII] line. The latter 
may require the selection of different spectral ranges depending of 
the source, because of the different location of the atmospheric ab¬ 


sorption features in different spectra (relative to the expected [CII] 
frequency). In the cases of BDF3299 and BDF521 the expected lo¬ 
cation of the [CII] line is in an excellent atmospheric region. In 
these cases we have measured the rms of the spectrum between - 
1000 km/s and -h300 km/s; outside this range the rms increases as 
a consequence of atmospheric absorption features. In the case of 
SDF46975, an atmospheric absorption feature is located just red- 
ward of the expected [CII] frequency from the Lya redshift; since 
we do not know exactly the expected location of the [CII] from the 
real accurate redshift, we conservatively estimate the noise from the 
spectral range corresponding to the atmospheric absorption feature. 

The absolute astrometric accuracy of the maps is given by the 
accuracy of the phase calibration and the accuracy of the absolute 
coordinates of the phase calibrator. The global absolute astrometric 
uncertainty is estimated to be about 0.15". In the case of BDF3299, 
two observations were performed with two different phase calibra¬ 
tors: in the two observations the continuum serendipitae have the 
same absolute position within less than 0.1", confirming the accu¬ 
racy of the absolute astrometric calibration of ALMA. 

The absolute astrometric accuracy is only relevant for the 
comparison with the Y-band images. In the case of BDF3299 the 
ALMA data reveal two serendipitous continuum sources detected 
at high level of significance ( |Carniani et al.|2015] l, both of which 
have a clear counterpart in the Y-band image or in the K-band im¬ 
age (which is cross-matched with the Y-band). However, both op¬ 
tical and near-IR counterparts are offset by about the same amount 
in the same direction relative to the ALMA continuum sources, and 
more specifically by ~ 0.4" towards the NW (PA~ 343°). Such a 
systemic offset is shown and discussed in AppendixSince these 
offsets are the same for both continuum sources we believe they 
are associated with an offset of the astrometry of the near-IR im¬ 
ages relative to the ALMA images (we shall mention that the ab¬ 
solute astrometry of the Y-band and K-band images is linked to the 
USNO-B catalog, which has an average absolute astrometric accu¬ 
racy of 0.25" at the epoch J2000, but variable across the sky with 
a dispersion about 0.12"). Therefore, when comparing with the Y- 
band image, we will apply to the Y-band image an offset of 0.4" to 
the SE (PA= 163°). 
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Figure 2. [CII]158jum luminosity versus SFR. The blue line shows the re¬ 
lation found by |De Looze et al.H2014) for local star forming galaxies and 
starbursts (non including ULIRGs), while the shaded region shows the Icr 
dispersion. The orange dashed line and hatched region show the relation 
and dispersion found by |De Looze et al.| (2014^ for local low metallicity 
dwarfs and irregulars. The black diamonds are galaxies at z ~ 5 - 6 from 
[Capak etn^(2015) and |Willott et al.|(2015b) and their average errorbar is 
shown in in the bottom-right of the diagram. The red symbols are the [CII] 
data for galaxies at z~7 (specifically at 6.8<z<ll): circles are our new ob¬ 
servations and squares are data from the literature (Ota et al.|2014||Schaerer| 
|et al.|2015bl|Gonzalez-L6pez et al.|201^ . 


3 [CII] EMISSION AT THE LOCATION OE THE 
OPTICAL SOURCES 

The [CII] line is not detected (at a confidence higher than 3cr) in 
any of the three sources at the location of their Y-band emission (i.e. 
within a radius 0.5" from the latter). We recall that in these sources 
the Y-band emission is partly contributed by UV stellar contin¬ 
uum emission and partly by Lycr emission (at least for BDF521 
and BDF3299, Lycr and UV continuum contribute to 1/3 and 2/3, 
respectively, of the observed Y-band fiux, |Vanzella et al.|2011| l. 

The spectra extracted from the central beam are shown in 
Fig.[^ along with the spectral distribution of the noise in the lower 
panels. 

To estimate the upper limit on the line flux we assume a line 
width of 100 km/s (we will see in the next section that this is the 
velocity width of the clump detected next to BDF3299, moreover 
this width is consistent with other lower redshift galaxies with low 
SFR, e.g. [Williams et al.|2014]|Capak et al.|2015| ). Note that other 
authors provide upper limits by assuming a line width of 40 km/s 
( jOta et al.|2014|[Gonzalez-L6pez et al.|2014| l; since these data will 
be used in our analysis, we have re-scaled their upper limits to a line 
width of 100 km/s for consistency (this is simply done by assuming 
that the noise on the line flux scales as the square root of the line 
width). The 3cr upper limits on the [CII] luminosity are given in 
Tab.[T] 

Since [CII] is primarily excited as a consequence of heating 
by soft-UV radiation (mostly in PDRs associated with star forming 
regions, but also in HIT regions or in the diffuse medium), many 


authors have identified a correlation between SFR and [CII] lumi¬ 
nosity (e.g jSargsyan et al.|2012[|De Looze et al.|2014j |. 

If the far-IR luminosity is used to measure the SFR, then the 
correlations seems to break down at high far-IR luminosities (in 
the ULIRGs regime) and it is not yet clear whether this is asso¬ 
ciated with an additional contribution to the far-IR from powerful 
AGNs, hence resulting in an inappropriate measurement of the star 
formation rate from the far-IR emission (Sargsyan et al.|201^|Mal-| 
jhotra et al.|2001j l, or associated to different physical conditions in 
ULIRGs, which may suppress the emission of [CII] (Kaufman et al.j 
|1999[ [Gracia-Carpio et al.||2011| l. However, our sources are cer¬ 
tainly not in the ULIRG luminosity range, so we will not consider 
the latter issue any further. 

Fig. shows the L([CII]) versus SFR relation for various 
classes of galaxies presented in jDe Looze et al.|(20I4j ). The blue 
line indicates the relation for local normal star forming and star- 
burst (non-ULIRG) galaxies, with the shaded blue region giving the 
+ lcr dispersion. The orange dashed line and hatched region show 
the relation and dispersion of local low metallicity galaxies (in the 
range ^ < Z/Z© < 1). For many of the low metallicity galaxies in 
the sample of jPe Looze et al.| ( |20I4j ) the L([CII])/SFR ratio tends 
to be even higher than in normal galaxies (this is mostly a conse¬ 
quence of [CII] being a coolant of the ISM), however at very low 
metallicities (Z < 0.1 Z©) the L([CII])/SFR start to decrease and 
becomes lower than in normal galaxies. 

Galaxies at 5<z<6.1, from [Capak et al.j (2015j ) and [Willott[ 
[et al.[ ( [^15b[ ), are shown with black diamonds and tend to follow 
the same [CII]/SFR relation as local galaxies. These galaxies, with 
SFR typically of a few/several times 10 M© yr"^ are not yet repre- 


sentative of the bulk of the population at these redshifts I 

[Salvaterra 

et al.[[2011|[Finkelstein et al.[[2012| [Dayal et al.[[20I4|| 

|Robertson 

et al.|2015>. 


The 3cr upper limits for “normal” galaxies at z~7 (those with 
SFR < 100 M©yr"^ i.e. not including the extreme QSO at z=7.08) 
are indicated with red symbols in Fig. Circles indicate our new 
observations, while squares are for additional galaxies at 6.8<z<ll 


from previously published studies ( 

Ota et al.|20I4| [Schaerer et al.[ 

[2015b|[Gonzalez-L6pez et al.|20I4 

(. At least three galaxies at z~7 


for which a meaningful constraint on L([CII]) has been obtained 
(including our ALMA deepest observation, BDF3299) have, at a 
given SFR, a 3cr upper limit on their [CII] luminosity significantly 
lower than observed in most galaxies at lower redshifts, even by 
including the low metallicity ones. 


However, we note that the lowest metallicity local galaxies 
in the [De Looze et al.] (2014[ ) sample (Z ~ 0.1 - 0.02 Z©) have 
[CII]/SFR values comparable with the 3cr upper limit on BDF3299 
(although they have much lower SFR). Therefore, in principle, such 
a [CII] line deficit in z~7 normal galaxies could be ascribed to low 
metallicity. Low metallicity of high-z galaxies has been a possible 
interpretation for the lack of [CII] detection in some star forming 
galaxies at z~6.5 (Ouchi et al.|2013[[Vallini et al.|2013[ l. Yet, one 
should bear in mind that the intensity of the [CII] line, as many 
other coolants of the ISM, does not scale linearly with metallicity 
(see e.g. Nagao et “^[20111 [Madden et al.[[20l H [Cormier et al.[ 


[2012[ [Rollig et al.[[2006[ l. Most importantly, we will show in the 

next section that [CII] emission is actually detected in a gas clump 
next to BDF3299. Therefore, while metallicity may play a role, it 
is probably not the only origin of the [CII] deficiency in the core of 
these primeval galaxies, as discussed in the following. 
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Figure 3. Top: Map of the [CII] emission obtained with a velocity interval 
of 100 km/s, centered at -64 km/s (relative to [CII] frequency expected 
from the Lya redshift, z=7.I09), in the ALMA data of BDF3299. Black 
contours are at levels of 2, 3 and 4 times the noise per beam in the same 
map, i.e. 6.4 mJy beam“^ km s“^. Grey contours are at the same levels, but 
for negative fluxes (note that -4cr is not reached anywhere). The ALMA 
beam is indicated in the bottom-right corner. The white cross indicates the 
location of the Y-band image centroid of the galaxy BDF3299. Label “A” 
indicates the contours associated with the [CII] detection, whose spectrum 
is shown in Fig.|^ while “Bl” and “B2” indicate two marginal detections 
discussed in the text. Bottom: zoom of the central 3" around the location of 
BDF3299. Note that the coordinates are relative to the phase center of the 
ALMA observation in 2013, which was centered on the nominal coordinates 
of the optical image (but see the text for a discussion on the astrometry of 
the original Y-band images and their shift relative to the ALMA astrometry, 
which has been corrected here). 
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Figure 4. Top: ALMA spectrum of the [CII] source (clump “A”), extracted 
from the an elliptical aperture obtained by fitting the emission in the [CII] 
map (the spectrum flux has been corrected for the flux lost outside the ex¬ 
traction aperture, see text). The gray shaded region indicates the part of 
the spectrum affected by higher noise because of atmospheric absorption. 
Bottom: Lycr profile from the optical spectrum obtained by |Vanzella et al.| 
j2011) . It should be noted the strongly asymmetric profile due to IGM ab¬ 
sorption of the blue side of the line. The green dotted line indicates the sky 
background spectrum. Velocities are relative to the peak of Lya (orange 
dot-dashed vertical line). The red dashed line indicates the velocity of the 
peak of the [CII] line. 


4 [CII] EMISSION NEXT TO BDF3299 AT Z=7.1 
4.1 Line detection 

In the ALMA cube of BDF3299, which is the deepest of our obser¬ 
vations, we have detected a line consistent with [CII] at the redshift, 
z=7.1, inferred from the Lycr of the galaxy (especially once IGM 
absorption of Lya is taken into account), but offset by 0.7" (i.e. 
4 kpc) relative to the location of the primary galaxy in the Y-band 
image (which traces UV-fLycr emission). 

The map of the line emission is shown in Fig.|^ extracted with 
a spectral width of 100 km/s and centered at -64 km/s, relative to 
the velocity scale set by the redshift inferred from Lycr. The white 
cross indicates the centroid of the Y-band image of BDF3299. 

The [CII] emitting clump is marked with an “A” and is located 
at RA(J2000)=22:28:12.325, DEC(J2000)=-35:10:00.64. Its emis¬ 
sion peak is detected at 4.5cr. However, the source is marginally 
resolved (along a Position Angle of -33.5°) hence the significance 
cannot be inferred simply by the peak surface density in units of 
fiux per beam. In this case, as for extended sources, to properly as¬ 
sess the significance of the detection, we have to extract the spec¬ 
trum of the source. We have first fitted a two-dimensional Gaus¬ 
sian to the [CII] emission, and then extracted the spectrum from an 
aperture corresponding to the section of the ellipsoid at half max¬ 
imum of its peak (i.e. an ellipses which has minor and major axes 
of 0.55" and 1.03" respectively). Since the extraction aperture has 
a size close to the beam, we have applied a flux aperture correction 






































6 Maiolino et al. 



1" = 5.2 kpc 

a 


1 

oO 



^ n 

0) o 

(0 
o 

Nd 







o 

o 

W 



1-2 


4 

1 B2 

Q 




§ 




CM 

“3 

> -4 

-M 

D 



0 

13 

Cd 

-5 

5 0 

» o 

-6 





5 4 3 2 

1 0 -1 

-2 


Relative J2000 Right Ascension (arcsec) 


1 

0.5" = 2.6 kpc 



0.5 

CO 

o 

4- 




c 

iS 0 

"d 

c 


jm 


"o 

(D 

Q 

J 

m 


o -0.5 

o 

o 

CM 



“D 

CD 

-1 

o 

0 




cn 

A 

O 

-1.5 





1.5 1 0.5 

o 

1 

p 

cn 

-1 


Relative J2000 Right Ascension (arcsec) 



Figure 5. Top panel: The color background image is the Y-band image of 
the BDF3299 field. The black cross indicates the location of the galaxy 
BDF3299 at z=7.109. We recall that at this redshift the Y-band samples both 
Lya and UV continuum. The green contours show the [CII] map with the 
same levels as (black contours) in Fig.|^ Label “A” indicates the contours 
associated with the [CII] detection, whose spectrum is shown in Fig. 
while “Bl” and “B2” indicate two marginal detections discussed in the text. 
The ALMA beam is shown in the bottom-right corner. Bottom panel: zoom 
of the central 3" around the location of BDF3299. Note that the coordinates 
are relative to the phase center of the ALMA observation in 2013, which 
was centered on the nominal coordinates of the optical image (but see the 
text for a discussion on the astrometry of the original Y-band images and 
their shift relative to the ALMA astrometry, which has been corrected here). 


(inferred from the flux extracted with the same aperture from the 
continuum serendipitae), which results to be a factor of 2.5. The re¬ 
sulting spectrum is shown in Fig. rebinned to 40 km/s, together 
with the optical spectrum of Lycr. The line has a peak at a veloc¬ 
ity of -71 km/s relative to the Lycr pealQ As already mentioned, 
such a velocity offset is expected. Indeed the Lycr is redshifted rel¬ 
ative to the systemic velocity as a consequence of IGM absorption 
(as clearly shown by its strongly asymmetric profile), and possibly 
also due to the effect of outflows (e.g. [Pettini et aL||2002| |, hence 
the real redshift of BDF3299 is certainly slightly lower than in¬ 
ferred from Lycr. Therefore, the blueshift by -71 km/s is fully con¬ 
sistent with the rest-frame of the galaxy. We have measured the 
significance of the line by measuring the noise in the spectrum be¬ 
tween -1000 km/s and -1-300 km/s; this is a region with uniform, 
low noise, while outside this spectral range atmospheric absorption 
makes the noise higher (see Fig[^. To avoid any potential baseline 
residual instability, we have also subtracted a continuum fitted on 
the same spectral window (which however is consistent with zero 
within the noise). By integrating the line under the shaded region 
in Fig. [^(-200 < V < 0 km/s), we obtain a significance of the 
detection of 7cr. A summary of the [CII] line properties of clump 
A is given in Tab. 

We note that the same map in Fig. [^reveals two potential ad¬ 
ditional marginal detections (Bl and B2), which are the only ones 
in the map with significance higher than 3cr within a radius of 6" 
and will be discussed later on. 

The spatial location of the [CII] detection relative to the Y- 
band source is shown in Fig. in which the green contours are 
the same as the (black) contours in Fig. while the colored back¬ 
ground image is the Y-band image. 

We note that at the location of the optical source the ALMA 
[CII] map shows a some signal at the 2cr level. However, this signal 
is very marginal and it will not be discussed any further. 


4.2 Reliability of the detection 

In this section we discuss more in detail the significance of the [CII] 
detection. 

First we have verified that the detection is not the result of a 
glitch in a single subset of the data. We have separated the visibil¬ 
ities obtained for BDF3299 in three data sets: the set of observa¬ 
tions obtained in November 2013, and we have divided into two 
halves the visibilities obtained in April 2014. These three data sub¬ 
sets have similar depths (though the data in 2013 are somewhat 
shallower than each of the two halves data taken in 2014). We have 
verified that in each of these three data sets a signal of the [CII] 
line is seen, at a flux level consistent with the value observed in 
the merged data, within the noise of each individual subset of data. 
More specifically, by only taking the (more sensitive) data taken in 
2014, the measured line flux is 46.0 + 7.7 mJy km s“^ If the 2014 
data (4 scans in total) are split in two groups of two scans each, then 
the line flux in each of the two scans is 44.7 + 14.0 mJy km s“^ and 
47.3 + 12.7 mJy km s"^ If only the (lower sensitivity) data taken in 
2013 are used, then the measured line flux is 49.9+17.9 mJy km s“^ 
(note that all of these fluxes, and associated errors, have been scaled 
by the same aperture correction factor of 2.5 discussed above). 

We have then checked whether negative sources are detected 
with the same significance or not. We have searched for additional 


^ This velocity is the average between the velocity estimated from a Gaus¬ 
sian fitting (-64 km/s) and the first moment of the line (-78 km/s). 
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Figure 6. Top: Distribution of positive (red) and negative (blue) detections 
in the SPWl of BDF3299, as a function of signal-to-noise ratio of the detec¬ 
tion, for lines integrated over the same velocity range as our [CII] detection, 
searched wihtin one primary beam and in the velocity range -1000 km/s < 
V < -1-300 km/s (i.e. within the range of good atmospheric transmission). 
Middle: cumulative distribution of positive (red) and negative (blue) de¬ 
tections, within the same ranges. Bottom: Distribution of positive (red) 
and negative (blue) detections within a radius of 3" from the position of 
BDF3299, and over the same velocity range -1000 km/s < v < -1-300 km/s. 


emission line detections within SPWl, within the spectral region 
between -1000 km/s and -1-300 km/s (i.e. the region in which the 
noise is low and not affected by atmospheric absorption features). 
We have adopted the same velocity width (100 km/s) as for the 
[CII] detection, to extract maps for the initial selection of putative 
sources. We have restricted the search to this criterion both because 
the sidelobes of the continuum serendipitae may affect different ve¬ 
locity ranges in a different way, and because we have verified that 
the noise does not scale exactly as (AV)"^/^, possibly because of 
some residual baseline problems (see also Vio et al., in prep, for 
a thorough analysis of the noise in the ALMA data). We have dis¬ 
carded putative detections (either positive or negative) for which: 
1) the signal is not seen at the expected level (within the noise) 
in the three subsets of data discussed above (i.e. the detection is 
mostly resulting from a glitch in one or two of the three subsets of 
data); 2) the measured size of the putative detection is smaller than 
the synthesized beam, or 3) larger than 1.5 times the beam (real 
sources are not expected to be very extended). The latter two re¬ 
quirements come from the fact that a source with size significantly 
smaller (or much larger) than the synthesized beam must be due 
to noise associated with individual antennae, or group of antennae 
(these sources of noise are not correlated and therefore should give 
features with angular profiles much different than the beam). Side- 
lobes of serendipitae, or from strong sources outside the field of 
view, would also introduce features with sizes different from the 
beam. However, to be on the safe side, we have also excluded the 
region within 2" from the 0.4 mJy continuum serendipitous source, 
since the uncleaned continuum map shows some low-level side- 
lobe residuals in this area. For each putative detection passing this 
screening we have extracted the spectra as for the [CII] detection 
and estimated the significance as discussed in Sect |4.1| (i.e. over the 
same velocity interval of 200 km/s as for our detection). The result¬ 
ing distribution of the positive (red) and negative (blue) detections 
is shown, as a function of the Signal-to-Noise ratio, in the top panel 
of Fig. The cumulative distribution of the positive and negative 
detections (i.e. the number of detection with significance above a 
given threshold) is shown in the middle panel of Fig.[^ 

It is interesting to note that the number of positive detections 
at > 3cr is 30% higher than the number of negative detections, sug¬ 
gesting that a significant fraction of the positive detections are real 
even down to the 3cr level. Focusing on the detections at > 5cr, we 
note that there are no negative detections, hence giving further con¬ 
fidence that our positive detection at 7cr next to the optical position 
of the galaxy (clump “A”), is real. 

If we restrict the search within a radius of 3" from the optical 
location of BDF3299 (but still ove the full velocity range -1000 < 
AV < -1-300 km s“^), then the resulting distribution of detections is 
shown in the bottom panel of Fig.[^ in which there are only three 
potential positive detections, and no negative detections, above 3cr. 
This further supports the reliability of our detection in clump “A”. 

4.3 Line identification 

The match of the line frequency with the one expected for [CII] at 
the redshift of BDF3299 (Fig. strongly supports this identifica¬ 
tion of the line. 

The identification with another transition associated with a 
foreground galaxy at lower redshift, which by chance happens to 
be at 0.7" from BDF3299, is extremely unlikely, since the sky den¬ 
sity of molecular and atomic lines in this band, and with fluxes 
comparable with the line detected in our data, is extremely low. 
Even in such unlikely case, the galaxy should be clearly seen in the 
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optical/near-IR images and in the thermal millimeter continuum. 
In the following we quantify these statements, hence excluding the 
line identification with a lower redshift transition. 

The CO rotational transitions are the brightest transitions 
that can be seen in low/intermediate redshift galaxies at this fre¬ 
quency. Other transitions are generally much fainter. [da Cunha et"^ 
( |2013bj ) have inferred the cumulative space density of CO transi¬ 
tions as a function of fiux detection threshold. 

We have considered their expectations including all CO ro¬ 
tational transitions from (7-6) to (2-1). jda Cunha et al.j ( |2013b| ) 
consider both the CO Spectral Line Energy Distribution (SLED) of 
the MW and the extreme CO SLED observed in the center of the 
starburst galaxy M82. We have first conservatively adopted the CO 
SLED of the M82 center; this is a conservative assumption, since 
the bulk of high-z galaxies have CO excitation closer to the MW 
(e.g. jPannerbauer et al.|2009| . With this assumption we obtain that 
the cumulative probability of finding a CO emitter within a radius 
of 1" of a given position, within a velocity range of +300 km s"^ 
from a given frequency in band 6, and with a flux brighter than 
30 mJy km/s (i.e. 50% lower than the fiux of the line detected by 
us), is about 10""^. 

If the CO SLED typical of the MW (which is closer to the CO 
SLED of the bulk of high-z galaxies) is taken, then the probability 
of chance CO detection drops to 2 10"^, within a velocity interval 
of +300 km s"^ and to 7 10"^, within a more reasonable velocity 
interval of +100 km s“^ 

In the calculations above we have conservatively assumed the 
putative foreground galaxies to have solar metallicity, which has 
been used to estimate the CO sky densities. However, high redshift 
galaxies are characterized by lower metallicity than local galaxies 
(e.g. [Maiolino et al.||2008] l. The CO-to-H 2 conversion factor de¬ 
pends steeply on metallicity ( [Bolatto et al.|2013) , with low metal¬ 
licity galaxies being much less efficient in emitting CO for a given 
molecular gas mass. Therefore, the probability of finding a fore¬ 
ground CO emitter in the vicinity of BDE3299, and at the same fre¬ 
quency expected for its [CII] emission, is even lower than inferred 
above. 

In addition, to such extremely low probabilities of having a 
chance detection of a foreground, low-z CO emitter, one should 
take into account that a putative CO emission, with the inten¬ 
sity observed by us, should be associated with a relatively bright 
galaxy. Indeed, according to the Schmidt-Kennicutt law, CO emis¬ 
sion tracing molecular gas in a galaxy, must be associated to a cer¬ 
tain SER. Assuming, for instance, that the observed line is CO(4- 
3) at z=l and conservatively assuming the S-K relation observed 
for local galaxies ( [Kennicutt & Evans|2012| ), one expects that the 
host galaxy should have a SER ^ 25Mo yr"k Such a SER should 
be associated with a galaxy with a brightness in the Y-band of 
Yab = 20 mag. However, there is not such bright counterpart at 
the location of the line detection. Beside BDE3299 itself, there is a 
faint source, marginally detected (Eig[^, but confirmed in new HST 
images (Castellano et al. in prep.), located 0.3" to the South-East, 
but which is seven magnitudes fainter than expected. The faintness 
in the Y band could in principle be associated with extreme dust 
extinction. However, a SER of 25Mo yr"^ should also result into a 
lOcr continuum detection in our ALMA data, while we do not even 
have a marginal detection at the location of the line. Similar or even 
more extreme expectations are obtained for other CO transitions or 
for any other submm molecular or atomic transition. 

These expected brightness in the Y-band and in the millime¬ 
ter should be considered as very conservative. Indeed, it has been 
shown that the S-K law evolves at high redshift, in the sense that 


for a given CO luminosity, or for a given gas content, high-z galax¬ 
ies are characterized by higher SER than local galaxies ( [Tacconi| 
[et al.[[20T3| [Santini et al.|[2014] ). High-z, low metallicity (faint) 
galaxies would be even more offset from the S-K relation, in the 
sense that their CO emission would be even fainter (because of 
the steep metallicity dependence of the CO-to-H 2 conversion fac¬ 
tor) or, equivalently, for a given CO luminosity the associated SER 
rate (hence Y-band or far-IR continuum luminosity) should be even 
higher than estimated above. 

Einally, the arguments above, based on the S-K law, apply to 
“normal” galaxies. Starburst galaxies have an efficiency of star for¬ 
mation (SER per unit gas mass) one order of magnitude higher 
( [Genzel et al.|[2010[ [Daddi et ^|2010j l. Therefore, depending on 
whether the putative foreground galaxy is “normal” or “starburst” 
the expected luminosity in any band should be a factor of a few, up 
to an order of magnitude, higher than estimated above. 

Summarizing, the identification of the detected line with an¬ 
other transition associated with a foreground galaxy is extremely 
remote, leaving [CII] at z=7.I07 as the only realistic viable inter¬ 
pretation. 

4.4 Interpretation of the [CII] emission 

[CII] emission offset relative to the UV and Lya emission was al¬ 
ready been seen in some other high redshift galaxies (e.g. [Williams I 
[et al.[[20T4| [Gallerani et al.|[20T^ , therefore our finding of off¬ 
set [CII] emission is not a novelty. However, previous cases were 
associated with peculiar objects, such as quasar host galaxies, or 
galaxies in an overdense environment and in the vicinity of a 
quasar. More recently such displacement between [CII] emisison 
and UV/Lyof emission has been observed also in high-z galaxies 
with more modest SER (a few/several times 10 M© yr"^), although 
still higher than the galaxies observed by us, and at lower redshifts, 
5<z<6.1 ( [Capak et al.|2"0T5l|Willott et al.|2015b| t. 

Interestingly, displaced [CII] emission in primeval galaxies, 
as observed in our ALMA data, is an expected consequence of 
the strong stellar feedback that, according to many models, should 
characterize galaxy formation in the early Universe (e.g. |Salvadori| 
|& Eerrara|2009[[Dayal et al.|2014[|Graziani et al.|2015| ». The [CII] 
line can be collisionally excited in many different ISM components 
(e.g. cold neutral medium, and high density photodissociation re¬ 
gions (PDRs) in molecular clouds). However, negative (thermal, 
radiative, mechanical) feedback from stars can efficiently disperse 
or disrupt nearby molecular clouds, therefore quenching their [CII] 
PDR emission. In this case, [CII] emission can only arise from ac¬ 
creting/satellite clumps of neutral gas displaced from the primary 
galaxy main body. 

To illustrate this scenario more quantitatively, i n Eig. [^ we 
show the model of primeval galaxies developed by [Valliniet al.| 
( [2013j ), but specifically tailored to a galaxy at z=7.1, with the same 
SER as BDE3299. In this model the gas in the main galaxy is com¬ 
pletely ionized, while molecular clouds (and the associated PDRs) 
are absent as a consequence of the strong stellar feedback (details 
on the model are summarized in Appendix [^. The black points 
show the distribution of young hot stars, i.e. what we observe as 
UV rest-frame continuum emission (we recall that the rest-frame 
UV-emission contributes to 2/3 of the observed Y-band emission). 
The color image shows the distribution of warm ionized gas, which 
can be observed as Lycr emission (note that the total observed Lycr 
emission can in principle be more extended as a consequence of 
scattering). The contours show the [CII]158yL/m emission associ¬ 
ated with the neutral gas (or mildly ionized gas). Clearly the central 
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Figure 7. a) Simulation of a primeval galaxy at z=7.1 with SFR similar to 
BDF3299. The distribution of column density of ionized gas (which emits 
Lya) is shown in colors (note that the observed diffuse Lycr emission can be 
more extended than the warm gas as a consequence of scattering). The black 
points show the distribution of young stars (UV continuum, which domi¬ 
nates the Y-band emission in our images). The contours show the emission 
of [CII] tracing the neutral gas (distributed in clumps orbiting/accreting the 
primary galaxy), b) Same as the top panel in which the image of the ionized 
gas has been convolved with the angular resolution of our Y-band images 
(0.5") and the [CII] map has been convolved with the ALMA beam of our 
data, c) Integrated [CII] spectrum from the whole simulation. 


galaxy is expected to emit strong UV continuum and Lycr, but being 
completely photoionized, it does not emit any [CII]. However, very 
interestingly, the model expects [CII] emission from gas clumps 
at a few/several kpc (~ 1" - 2") from the primary galaxy. These 
are satellite gas clumps in the process of accreting onto the pri¬ 
mary galaxy, which have survived photoionization, owing to their 
distance from the source of UV photons (and because a significant 
fraction of UV photons are absorbed internally by the ISM of the 
primary galaxy). 

Such [CII] emission in the satellite, accreting gas clumps 
is expected to be faint, but our deepest observation, the one of 
BDF3299, does have the sensitivity to detect it. 

In Fig.[^ we have smoothed the simulation to the angular res¬ 
olution of our optical observations (0.5" seeing) and to the ALMA 
beam of our observations. The simulation is not meant, by any 
means, to exactly reproduce our observations (this would require 
a huge number of simulated objects with the same detailed mod¬ 
elling of the photoionization and radiative transfer in the ISM and 
IGM shown here, and then finding the best matching case). How¬ 
ever, Fig. illustrates that the offset between Y-band emission 
and [CII] emission is indeed expected to be resolved for various 
[CII] clumps (and especially for the brightest one) with the angu¬ 
lar resolution delivered by our ALMA observation. Lower angular 
resolution observations would probably hamper the capability of 
resolving the spatial offset and would probably associate the [CII] 
emission with the primary UV-Lycr emitting galaxy. 

The simulation also highlights that more than one [CII] clump 
may be present around BDF3299, out to 3" and possibly even at 
larger distances (3" is the maximum distance probed by the simu¬ 
lation box, so potential clumps at larger distances escape the sim¬ 
ulation). If present, these clumps are probably below our detection 
limit and may also suffer beam dilution (many of the isolated small 
clumps disappear in the smoothed simulated map). Actually, in the 
ALMA cube of BDF3299 we find other positive [CII] detections, 
at lower significance than clump “A”, but in excess relative to the 
“negative” detections, hence suggesting that they are not spurious 
detections. Two of such marginal detections are marked as B1 and 
B2 in Fig.[^ These are the only two features with a peak flux de¬ 
tected at a significance higher than 3cr within a radius of 6" (half 
a beam radius) from the primary galaxy, in the same map, centered 
at -64 km s"^ Certainly, deeper and/or higher angular resolution 
observations are required to further confirm these marginal detec¬ 
tions. 

In principle the line width and spatial extension can pro¬ 
vide information on the dynamical mass of the blob. Determin¬ 
ing a dynamical mass of the [CII] clump detected in the vicinity 
of BDF3299 is not easy, both because we do not know whether 
the line width is tracing a rotating system or a velocity dispersion 
dominated system, and because the [CII] map is only marginally 
resolved, hence the intrinsic dimension (beam deconvolved) is dif¬ 
ficult to estimate. If we assume that the system is rotationally sup¬ 
ported and with a radius of about 1.5 kpc (inferred from the Gaus¬ 
sian fit deconvolved from the beam in quadrature), we obtain a 
rough estimate of the clump dynamical mass of about 5 x 10^ Mq. 
If the system is not virialised, then this is actually an upper limit 
to the dynamical mass. Although all the dynamical mass estimate 
is extremely uncertain, it suggests that indeed the [CII] clump is a 
small satellite of the primary galaxy. 

According to the model, in these satellite/accreting clumps the 
[CII] excitation is not expected to primarily originate from in-situ 
star formation, but mostly as a consequence of heating from the 
soft-UV irradiation from the central primary galaxy. This is con- 
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sistent with the fact that we do not detect a clear Y-band source 
(i.e. Lycr+UV continuum) at the location of the [CII] clump. The 
vicinity of BDF3299 (Yab = 26.15) at 0.7", and of another faint 
foreground galaxy at 0.4" (Yab ~ 27.5, see next section), prevent 
us to set tight upper limits on the Y-band emission associated with 
the [CII] clump, but it is certainly at least 3-4 times fainter than the 
primary galaxy. According to the [CIIJ-SFR relation for local (and 
intermediate redshift) galaxies ( |De Looze et al.||2014| ), the [CII] 
luminosity of the clump would imply a SFR ~ 6 M© yr“\ i.e. com¬ 
parable to the main galaxy BDF3299, i.e. the [CII] clump should 
have similar or higher Y-band flux. The inconsistency between the 
[CII] clump and the local [CII]-SFR relation is also illustrated in 
Fig. which also shows the location of the [CII] clump, unde¬ 
tected in Y-band, on the [CII]-SFR diagram. The clump is clearly 
inconsistent with the distribution of local galaxies, but it could still 
be marginally consistent with the expectations for low metallicity 
galaxies. Therefore, we cannot exclude some contribution to the 
[CII] excitation from in-situ star formation, but we can also state 
that the non-detection of a clear Y-band counterpart is certainly 
consistent with the models expectations. 

An alternative interpretation could be that the rest-frame UV 
emission associated with putative star formation in the [CII] clump 
is signiflcantly obscured by dust. However, in this case one would 
expect associated dust emission at the level of ~ 15 jiSy, which is 
~ 2 times higher than the rms of the continuum map, while we 
currently do not even have a mariginal continuum detection. We 
also note that the [CII] clump is not detected at longer wavelengths 
(Jab > 26.5, Kab > 26.0). Signiflcant dust content responsible for 
signiflcant extinction has also been generally excluded by various 
works on such high redshift galaxies {[W alter et al.|2012|[Ota et al.| 

|2014[ ISchaerer et al.|2015bl |Capak et al.|2015| l, as also discussed 

in this paper (in Sect. [^. However, the current data are not deep 
enough to completely exclude the scenario of dust obscuration. Ad¬ 
ditional data are certainly required to further dismiss or validate this 
possibility. Deeper HST data are currently being obtained and fu¬ 
ture JWST data will certainly provide a deflnitive answer. 

Obviously the clumps need some minimum level of metal en¬ 
richment to emit [CII]. The model suggests a metallicity of about 
0.05 Zq. However, the metallicity is very difficult to constrain by 
using a single emission line and especially based on [CII], since 
this one of the primary coolants of the ISM. Therefore, the metal¬ 
licity of the neutral gas clumps is highly uncertain and additional 
data are required to properly constrain it. Here we only note that 
a small level of metal enrichment in this clump may result either 
from pre-enrichment as a consequence of metals expelled from the 
primary galaxy, or from in-situ production. Indeed, even a very low 
SFR of only 0.2 Mo/yr, certainly undetectable by our observations 
(as discussed above), would be enough to enrich this small clump 
to a metallicity of ~0.05 Z© in less than 100 Myr. 

Fig. shows the integrated [CII] spectrum from the simu¬ 
lation. We note that the integrated [CII] spectrum inferred from 
the simulation, shown in Fig. [^, expects additional narrow emis¬ 
sion [CII] lines emitted from small accreting/satellite clumps il¬ 
luminated by the primary source. We have tentatively identified 
narrower (10-20 km/s) [CII] emission sources in the vicinity of 
BDF3299, but which need to be confirmed with deeper and higher 
angular resolution observations. 

The ALMA data of the other two sources in our sample 
(BDF521 and SDF46975) do not reveal clear evidence of [CII] 
emission in the vicinity of the UV-Lycr primary galaxy. However, 
the latter data are significantly shallower than the observation of 
BDF3299 (at least a factor of two in sensitivity). The same line 


seen in BDF3299 would be undetectable in the lower S/N observa¬ 
tions of BDF521 and SDF46975. 

Moreover, these two sources are brighter in the UV relative to 
BDF3299, hence their photoionization effect may extend to larger 
distances, therefore preventing neutral clouds to survive in the 
proximity of the galaxy. 


4.5 Possible lensing from foreground galaxies 


We already mentioned that a faint galaxy, marginally detected in 
the Y-band (but confirmed by recent HST data, Castellano et al. in 
prep.), is found at 0.3" to the SE of the [CII] emitting clump. The 
galaxy is marginally detected also in the V band, so it must be a 
foreground galaxy, certainly not associated with the z=7.1 system. 
Proper constraints on the redshift of this foreground galaxy will 
come from additional HST observations, as well as from forthcom¬ 
ing KMOS-VLT observations. Here we note that finding a fore¬ 
ground galaxy so close to the ALMA line detection is intriguing. 

The possibility that the line is actually CO at low redshift, 
and associated with this foreground galaxy, was already discussed 
and discarded above: both the extremely low sky density of CO 
emitters and the extremely faint counterpart (both in the Y-band and 
in the millimeter continuum) make this scenario implausible. Here 
we only add, to the arguments discussed above, that the putative 
CO emission would be offset from the Y-band counterpart, which 
would imply molecular gas offset by several kpc from the galaxy 
center, which would be at odds with any known galaxy, where the 
molecular gas is found within the central few kpc of galaxies. 

It is possible that the proximity of this source is simply re¬ 
sulting from a chance superposition. The source has a magnitude 
AB ~ 27 and den sity of sources at such faint fluxes is about 
100 arcmin"^mag"^ (Guo et al. 2013| ). The probability of randomly 
finding a galaxy with this magnitude, or brighter, within a radius of 
0.3" of a given position is less than 10“^. This probability is larger 
than that estimated for the chance CO detection discussed above, 
but still relatively small. 

Another possible scenario is that the foreground low-z galaxy 
is gravitationally magnifying the [CII] emission at z=7.1. Indeed, 
according to models, the primary galaxy at z=7.1 should be sur¬ 
rounded by several [CII] emitting clumps. Although these clumps 
are abundant within a radius of ~ 6" of the primary galaxy (30 kpc), 
the [CII] emission is probably below the detection threshold of our 
observations for most of them, especially as a consequence of their 
low metallicity. If a foreground galaxy happens to be in the held, 
this can introduce an even mild gravitational lensing that can boost 
the [CII] flux above our detection threshold. 

This scenario would naturally explain the proximity of [CII] 
emission at z=7.1 with a foreground galaxy, without invoking very 
unlikely chance overpositions. Indeed, in this scenario the several 
[CII] clumps become visible (their emission is boosted above the 
detection threshold) only when a foreground galaxy is located next 
to them in projection, so that some gravitational lensing can occur. 

The effect is particularly effective if the foreground galaxy is 
at z=1.5, which is the redshift maximizing the lensing magnifica¬ 
tion of radiation at z=7.1. Unfortunately we do not have elements 
to test this scenario more quantitatively, since we do not yet have 
information on the redshift, mass and light profile of the putative 
lensing galaxy and we do not have high enough resolution ALMA 
observations to properly investigate the morphology of the putative 
lensed emission. Therefore, for the time being this has to be re¬ 
garded as a speculation, which would explain the observations, but 
to be tested with additional data. 
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Figure 8. Far-IR (/Irest = 160//m) to far-UV (/Irest = 1500A) luminosity 
ratio for local galaxies as a function of metallicity (black symbols). The 
blue dashed line is a linear fit to the data. The colored solid lines indicate 
the upper limits on (yFy)FiR/(yFv)uv inferred for the sample of our galax¬ 
ies at z~7, specifically: red for BDF3299, orange for SDF46975, green for 
BDF521. 


However, in support of this lensing scenario, we note that the 
other two marginal [CII] detections in the vicinity of BDF3299 (B1 
and B2 in Fig[^, are also found within 0.4" of two foreground 
galaxies, as seen in Fig[^ The galaxy next to B1 is only marginally 
detected in our Y-band image, but it is confirmed in new HST data 
(Castellano et al. in prep.). The galaxy next to B2 is clearly detected 
and the data currently available give an approximate photometric 
redshift Zphot ~ 1, which is consistent with the lensing scenario. 

The cumulative probability of having, by random chance su¬ 
perposition, three Y-band detections within 0.4" of three given po¬ 
sitions (our [CII] detection “A” and the other two marginal [CII] 
detections “Bl” and “B2”) is 10"^. The lensing scenario would in¬ 
stead naturally explain the co-location of [CII] emission and fore¬ 
ground galaxies, without having to invoke extremely low chance 
probabilities. 

Certainly additional data are required to confirm and quantita¬ 
tively test the lensing scenario. 


5 CONTINUUM DUST EMISSION 

While we have found a few serendipitous continuum detections 
( [Carniani et al.|2015] l, the continuum dust thermal emission is not 
detected in any of our galaxies. This result is in line with previous 
results at high redshift ( [Walter et al.|2012[[Ota et al.|2014[|Schaerer| 
jet al.|2015b] ), which have found that distant “normal” star forming 
galaxies, close to the re-ionization epoch, have a rest-frame far- 
IR/UV luminosity ratio significantly lower relative to the bulk of 
local star forming galaxies, with upper limits on the far-IR emis¬ 
sion typically consistent with local low-metallicity dwarf or irreg¬ 
ular galaxies. This suggests that these early systems are character¬ 
ized by a small content of dust (which should reprocess UV light 
into IR thermal emission), both because of their low metallicity and 
because of shortage of time to produce dust (jValiante et al.|200^ 
[2011 [[Schneider et al.||2014[|Hirashita et al.||20141|Nozawa et al.j 


[2007[ ). Our non detections go in the same direction, but achieving 
even tighter upper limits on the far-IR/UV luminosity ratio. 

Fig. shows the ratio (yFy)FiR/(yFy)uv versus galaxy metal¬ 
licity for local galaxies, where (yFy)uv is estimated at Tj-est =1530A 
(i.e. GALEX FUV band, for galaxies in the local Universe), and 
(yFy)FiR is estimated at Trest =l60fim (i.e. Herschel-PACS longest 
wavelength band, for galaxies in the local Universe). Data are taken 
from [Dale et al.[ ( |?007[ ), [Dale et al.[ ( [T012[ |, [Ota et al.[ ( |20T4^ and 
priv. comm.) and from [Moustakas et al.[([2010[l (in the latter we 
have used the calibrations from jPilyugin & Thuan[2005[ ). We have 
applied the correction to the far-IR fluxes, as discussed in [Ota et ^ 
( [20I4[ ) and [da Cunha et al.[ ( [2013a[ ), to take into account the fact 
that at high redshift the detectability of thermal emission is more 
difficult because of the CMB temperature being non-negligible rel¬ 
ative to the dust temperature in galaxies (for the sake of simplicity 
we have adopted the same correction assuming all three sources at 
z-1). For local galaxies, there is a clear trend of (yFy)FiR/(yFy)uv 
to increase with metallicity, implying that a larger fraction of the 
UV light is reprocessed by dust at high metallicities. This is very 
likely a consequence of the dust content scaling with metallicity 
(e.g. [Remy-Ruyer et al.[2014[[Draine et al.[2007[ ). The blue dashed 
line is a linear fit to the data in the log-log plane. The 3cr upper 
limits on the (yFy)FiR/(yFy)uv on z~7 sources are indicated with 
horizonal colored lines (see legend). Clearly, galaxies at z~7 are 
inconsistent with the UV-to-IR SED of typical local galaxies, with 
solar-like metallicities, but are consistent with metal-poor (hence 
dust-poor) local dwarf galaxies, with metallicities about 7 times 
lower than solar. Yet, while this suggests that the z~7 galaxies in 
our sample must have metallicities similar or lower than local metal 
poor dwarfs, the latter typically have SER’s one or two orders of 
magnitude lower than the sources observed by us. Therefore, the 
galaxies observed by us at z=7.1 are probably scaled up versions 
(by at least an order of magnitude) of local low metallicity dwarfs. 

Constraints on the dust mass are in principle physically more 
interesting. However, it is not possible to obtain an upper limit on 
the dust mass based on a single upper limit on the rest-frame far-IR 
or submm emission, since the dust temperature is not known. [Ota[ 
[et al.[ ( [20T4| ) assumes a dust temperature of 27.6 K (and an emissiv- 
ity index = 1.5) resulting from the average of 17 local dwarf and 
irregular galaxies. [Schaerer et al.[ ( [2015b] ) assume a range of tem¬ 
peratures between 25 K and 45 K. We follow a similar approach as 
these papers to infer an upper limit on the dust mass from our 3cr 
upper limits at dobs = 1.2 mm, corrected for the effects of CMB. 
Depending on the assumed dust temperatures we obtain upper lim¬ 
its on the dust mass in the range between a few times 10^ M© and 
several times 10^ M©, as detailed in Tab.|^ We note that such low 
dust content is in agreement with the finding that the SER inferred 
from Lycr is consistent, within uncertainties, with the SER inferred 
from the UV emission ( [Vanzella et al.[20TT] |. 

In the early universe the dust content is not only determined 
by the metallicity, but also by timescale issues associated with the 
dust production mechanisms. Indeed, AGB stars, which are thought 
to be the main factory of dust in the loca l universe (e.g. [Ventura[ 
[et al.[[2012b[a[ [Pi Criscienzo et al.[[2013[ [Schneider et al. [[201 41 17 

require a timescale of several 100 Myr to evolve and to contribute 
significantly to the dust content, hence a timescale comparable to 
the age of the universe at z~7. Dust forming in the ejecta of core¬ 
collapse supernovae can be an additional source of dust on very 
short timescales, and indications that SN dust can indeed contribute 
significantly to the overall dust budget in the early universe, has 
been found by observational studies (e. g. [Maiolino et aT][2004[ 
[Gallerani et al.|2010| [Stratta et al.|2007[ [2011[ |. Extensive models 
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have been developed to investigate the expected dust evolution in 
primordial galaxies (e.g. Valiante et al.|2009]|2011[[Hirashita et al.| 


|2014[[Nozawa et al.|2015[|Gall et al.|2011 1. Our tight upper limits 

on the dust mass in these early systems further suggest that at these 
early epochs primeval galaxies have little time to produce masses 
of dust comparable with lower redshift galaxies. A more quantita¬ 
tive comparison with models of dust formation in the early universe 
goes beyond the scope of this paper and will be discussed in other 
forthcoming papers (Valiante et al. in prep., Ferrara et al. in prep.). 


6 CONCLUSIONS 

We have presented new ALMA millimeter observations aimed at 
constraining the [CII]158yum and thermal dust continuum emission 
in three primeval galaxies at z~7, whose redshift had been spectro¬ 
scopically confirmed through the detection of (asymmetric) Lycr. 
In contrast to many previous millimeter observations of high red¬ 
shift galaxies, which have targeted galaxies with extreme star for¬ 
mation rates (100-1000 M© yr"^), the galaxies in our sample have 
SFR ~ 5 - 15Mo yr"^ more typical of the bulk of galaxies at these 
early epochs. 

At the location of the optical (Y-band) counterpart, which 
samples the rest-frame Lycr and UV emission, no [CII]158yum line 
emission is detected. For some of the galaxies the 3cr upper limit 
on the [CII] luminosity is significantly lower than what is observed 
in local galaxies with similar star formation rates, even including 
many metal poor local galaxies The upper limit on the [CII]/SFR 
ratio in these galaxies is also lower than observed in galaxies at 
z ~ 5 - 6. 

However, in the deepest of our ALMA observations, targeting 
the highest redshift galaxy of our sample (BDF3299 at z=7.109) we 
do detect [CII] emission, slightly spatially offset, by 0.7" (4 kpc), 
relative to the optical (Y-band) counterpart. 

The line is detected at 7cr. The level of confidence is verified 
by checking the detection in subsets of the same data, by checking 
the absence of any negative detection at confidence level higher 
than 5cr in datacube over the whole ALMA primary beam and over 
a broad frequency interval, as well as the absence of any negative 
detection at a confidence level higher than 3cr within a radius of 3" 
from the optical source. 

The detection of [CII] spatially offset from the optical galaxy, 
along with the absence of [CII] emission at the location of the op¬ 
tical counterpart, is in nice agreement with the predictions of re¬ 
cent models of galaxy formation in the early universe. According 
to these models, strong negative stellar feedback in these primeval 
systems quickly destroys or disperses molecular clouds, hence 
quenching the emission of [CII] (the bulk of carbon is therefore in 
higher ionization states). However, [CII] emission arises from ac¬ 
creting/satellite clumps of neutral gas displaced from the primary 
galaxy main body. 

Therefore, our ALMA observations are directly probing such 
early phases of galaxy formation, in which both stellar feedback 
and gas accretion are at work. 

In the accreting clump [CII] emission is likely resulting pri¬ 
marily as a consequence of heating by the strong UV radiation from 
the primary galaxy. However, some contribution to the [CII] exci¬ 
tation from in-situ star formation cannot be excluded. 

At larger distances from the same galaxy (~ 2" - 4", i.e. ~ 
10-20 kpc) we obtain two additional marginal [CII] detections at 
the redshift consistent with Lycr. These marginal detections need to 
be confirmed with additional data, but their presence would be in 


line with the expectations of the models discussed above, according 
to which several satellites clumps of neutral gas should be emitting 
faint [CII] line in the surrounding of the primary primeval galaxy. 

It is interesting that both the [CII] detection and the two 
marginal detections are located in the vicinity (within 0.3" - 0.5") 
of foreground galaxies at lower redshifts. We show that the asso¬ 
ciation of the millimeter lines with other transitions at lower red¬ 
shift (in particular CO transitions) is extremely remote and incon¬ 
sistent with the properties of the counterparts (in terms of SFR). 
Chance superposition coincidence is possible, but very unlikely. We 
suggest, as an alternative explanation, that gravitational lensing by 
foreground galaxies may boost the fiux of the [CII] clumps located 
in the vicinity of the primary z=7.1 galaxy, hence bringing them 
above the detection threshold. Testing this scenario more quanti¬ 
tatively requires a deeper knowledge of the properties of the fore¬ 
ground galaxies (redshifts, masses, mass profile distribution, etc...) 
and higher angular resolution data of the [CII] emission. However, 
if confirmed, this scenario would indicate that lensing from fore¬ 
ground galaxies at z~l-2 could be an effective way to detect faint 
neutral clumps of accreting gas in the vicinity of primeval galaxies 
around re-ionization. 

Continuum thermal emission is not detected in any of the three 
galaxies of our sample. The 3cr upper limit on the (rest-frame) IR- 
to-UV continuum emission ratio is consistent with the value ob¬ 
served for the lowest metallicity local galaxies, in which the dust 
content is very low. Constraints on dust masses are difficult to infer, 
because only one single photometric upper limit per source is avail¬ 
able for the IR SED, hence no constraints on the dust temperature 
are available. By assuming the dust temperature typical of local low 
metallicity galaxies, we obtain 3cr upper limits on the dust masses 
of a few times lO^Mdust- The latter confirms that these early systems 
have a very low dust content, both because of low metallicity and 
because of shortage of time to produce dust at these early epochs 
through some of the standard dust production channels. 
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Table 2. Parameters of the [CII] detection next to BDF3299 in clump “A”. 


Vobs([CII])^ 

[GHz] 

Z[CII] 

AV'’ 

[km s“*] 

FWHM^ 
[km s“^] 

F([CII]) 
[mJy km s“^] 

L([CII]) 
[10’ Lo] 

234.43 

7.107 

-71+10 

102+21 

48.6+6.9 

5.9+0.8 


NotesObserved central frequency in LSRK. ^ Velocity offset relative to the peak of Lycr (AV) and full width half maximum 
(FWHM) are calculated both through the moment analysis and through a Gaussian fit; the values reported in the table are 
average values. 


Table 3. Constraints on the continuum upper limits. 


Name 

(vFy)FiR/(yFy)uv" 

Mdust(T = 27.6K) 
[lO’Mo] 

Mdust(T = 45K) 
[lO’Mo] 

BDF-3299 

<0.21 

<2.1 

<0.32 

BDF-512 

<0.44 

<5.2 

<0.77 

SDF-46975 

<0.23 

<5.8 

<0.87 


Notes:^ (vFy)FiR is estimated at drest = I60jum, (yFv)uv is estimated at drest = 1500A. 
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APPENDIX A: CONTINUUM SERENDIPITAE AND 
REGISTRATION WITH THE OPTICAL IMAGES 

The continuum map of BDP3299 reveals three serendipitous de¬ 
tections ( [Camiani et al.|2015] ): two detections within the primary 
beam (one at ~ 40cr and another one at ~ 6cr) and a third deetction 
outside the primary beam (at ~ 4cr). Pig. |A1| shows the continuum 
maps of the two detections within the primary beam. Both detec¬ 
tions have counterpart in the near-IR images (either Y-band or K- 
band). However, the centroid of the near-IR counterparts is slightly 
offset relative to the millimeter emission in the ALMA map. This is 
visible in Pig. |A1 [ where the white crosses show the location of the 
near-IR counterparts (in images which are not yet registered as dis¬ 
cussed in sect[^. Clearly both counterparts are offset by a similar 
amount and in the same direction and more specifically by about 
0.4" towards the NW (PA~ 343°). Since these offsets are the same 
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Figure Al. Top: Continuum map of the brigest serendipitous source de¬ 
tected in the field of BDF3299. Contours are at of lOcr, 20cr, 30cr, 40cr, 
where the rms=7.8//Jy. Bottom: Continuum map of the second brightest 
serendipitous source in the same field. Contours are at of 2cr, 4cr, 6cr. In 
both panels the white cross indicates the location of the near-IR counter¬ 
part. 


for both continuum sources we believe they are associated with an 
offset of the astrometry of the near-IR images relative to the ALMA 
images. We have therefore corrected for such an offset when com¬ 
paring the ALMA and near-IR images. 


APPENDIX B: THE THEORETICAL MODEL 

A detailed description of the model and simulations used in this 
paper is given in [Vallini et al.j p013[ ). In this appendix we only 
summarize some of the key points. 
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The |Vallini et al.|j2013^ model is based on cosmological SPH 
hydrodynamic simulations performed with GADGET-2 ISpringell 


comoving volume with 2 X 512^ baryonic-hdark matter particles, 
giving a mass resolution of 1.32 X 10^ M© for baryons and 6.68 x 
10^ Mq for dark matter. We select a snapshot at redshift z~7, and 
we identify the most massive halo (total mass Mh = 1.17x10^^ M©, 
rvir ~ 20 kpc) by using a Friend-of-Friend algorithm. We then 
zoom the simulation by selecting a (0.65 h^Mpc)^ comoving vol¬ 
ume around the center of the halo, and post-process UV radiative 
transfer (RT) using LICORICE ( |Baek et al.|2009| l. Gas properties 
are resolved on a fixed grid with a resolution of 60 pc. To define 
the position of the ionizing sources we assume that stars form in 
those cells characterized by a gas density p > 1 cm"^ in order to 
reproduce the typical size (~l-2 kpc) of star forming regions at z 
^ 6-7 pouwens et al.|2004||0uchi et al.|2009) , as inferred by UV 
continuum images. The population synthesis code STARBURST99 
peitherer et al.|199^ is used to obtain the ionizing spectrum of the 
galaxy. 

The simulation is complemented with a sub-grid model tak¬ 
ing into account the cooling and heating processes producing two- 
phase thermal structure of the neutral gas in the ISM (as in |Wolfire| 
|et al.|pV^ |2003| ). The density (ncNM, hwnm) and temperature 
(Tcnm, Twnm) of the the cold (CNM) and warm (WNM) neutral 
phase in each cell of the simulation is calculated as a function of : 

1) the gas metallicity, Z, determining the coolants abundance, and 

2) the Far- Ultraviolet (FUV) flux. Go, in the Habing (6-13.6 eV) 
band, controlling the photoelectric heating produced by dust grains. 
Go scales with the SFR and depends on the assumed age of the stel¬ 
lar population. The density and temperature of the CNM and WNM 
allow us to compute [C II] emissivity from each cell of the simu¬ 
lated galaxy as in Eq. (2) in |Vallini et al.| ( |201^ . 

Variations of such model in which molecular gas and PDR 
are present in the main galaxy are presented in (Vallini et al. 2015, 
submitted). 
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(20051. The original simulation reproduces a (10 h^Mpc)^ 























